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Abstract 



We consider the thermally activated fragmentation of a homopolymer chain. In our simple 
' model the dynamics of the intact chain is a Rouse one until a bond breaks and bond breakdown is 

oo : 

considered as a first passage problem over a barrier to an absorbing boundary. Using the framework 
of the Wilemski-Fixman approximation we calculate activation times of individual bonds for free 
§ | and grafted chains. We show that these times crucially depend on the length of the chain and 

the location of the bond yielding a minimum at the free chain ends. Theoretical findings are 

0" 

i i qualitatively confirmed by Brownian dynamics simulations. 
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The problem of thermally activated chain fragmentation (thermolysis) is of fundamen- 
tal interest in understanding the degradation and stabilization properties of polymers 
Whenever the corresponding fragmentation kernel (i.e. the probability per unit time for a 
bond at a given position to break) is known as the function of the breakdown position, time, 
chain's length etc., the overall fragmentation process is well-described, and the distribution 
of fragments at whatever time can be obtained by the solution of the corresponding kinetic 
equation [2]. Here the fragmentation kernel is the input into the universal theory, and many 
model forms for such kernels were formulated on the basis of parameterizing experimental 
observations or as simple analytical examples. However, up to our notion, the question of 
how does the corresponding kernel follows from the single polymer chain's dynamics was 
hardly considered. The only example known to us is contained in [2I which however deals 
with a model whose relation to the standard polymer dynamics ones is not quite evident. 

In what follows we discuss the thermally activated breakdown within a model which 
assumes that the dynamics of the chain is a Rouse one j^, Q (as long as it doesn't break), 
i.e. we disregard excluded volume effects, which is a reasonable assumption in the case of 
tagged chains in melts and concentrated solutions, as well as hydrodynamical interactions. 
We moreover assume that that the breakdown of the bond (represented as the breakable 
harmonic spring) takes place as soon as its elongation achieves the preassigned value (equal 
for all bonds in the chain). The possibility of reestablishing the bond after breakdown (defect 
healing) is neglected. We consider the relevant situations of free and grafted chains. 

The thermally activated bond breakdown is essentially an example of an intrachain chem- 
ical reaction. However, contrary to e.g. polymer cyclization B, 0, y, y, 0, y, q which 

was considered in quite a detail, this one was hardly tackled. The assumption that the bond 
breaks when achieving the given elongation simplifies the description, and corresponds to 
assuming the corresponding reaction to be purely diffusion-controlled reaction on a contact. 
Such problems can be casted into mathematical form of the first passage problem over a 
barrier to an absorbing boundary. Although the formulation of the overall problem is ex- 
tremely simple, its solution is not, since the projection of the overall chain motion onto the 



reaction coordinate makes the corresponding diffusion strongly non-Markovian [13j. This 
strongly non-Markovian nature of the problem involving multiple characteristic timescales 
reflects the fact that the reaction essentially takes place in a many-particle system. 

In recent years theoretical methods have emerged to treat diffusion-controlled reactions 
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among sites attached to polymers. Pioneering works go back to Wilemski and Fixman 
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15] and conceptual advances were made by Doi Q], de Gennes [3] and others 13]. But 
for all that, except for some special cases [7j, the analytical theories of reaction diffusion in 
polymer physics fail to give an exact description of the reaction rates and rely on additional 
assumptions [llj. However, as we proceed to show, the outcome of theoretical considerations 
within the framework of Wilemski and Fixman agrees qualitatively very good with the results 
obtained in Brownia n dy namics simulations. Distinct from numerous studies on end-chain 
reactions 111 . Il2| or the studies on interior loop forming reactions [l7j, we 

focus on the related but somewhat different problem of thermal activation of bonds, i.e., the 
first passage problem of nearest neighbor monomer distances. 

As we proceed to show the dissociation dynamics of the chain strongly depends on the 
location of the bond within the chain and the size of the system. Although the equilibrium 
distributions and activation barriers are the same for all the bonds, their activation times are 
not. At the free ends of the chain the first passage times are substantially lower compared 
to bonds in the middle of the chain (for a free chain) or at its grafted end (for a grafted 
one). Thus, the thermally activated fragmentation is expected to happen mostly at the 
chain ends. A similar behavior was found experimentally [li|. Furthermore it was shown 
that the forced rupture of adhesive contacts is strongly influenced by chain dynamics for 



undercritical forcing 



ISO. 



The paper is structured as follows: In the next section we introduce our model and discuss 
its breakdown properties within the simple one- dimensional setup. We recall the dynamics 
of the underlying polymer model and in Sec. [Ill an approach to calculate first passage times 
in diffusion controlled reactions based on Wilemski-Fixman approximation. In Sec. II III we 
study the impact of the location of the bond on its activation time for various chain lengths 
and activation energies. In Sec. II VI we generalize our approach to the experimentally relevant 
three-dimensional chain. Finally we summarize our results. 



3 



I. THE MODEL 



We first consider a one dimensional chain of N + 1 monomers with coordinates go, Qn- 
The interaction potential is given by 



U(q ,...,q N ) = -K^2 Rijqiqj 



(1) 



'■j 



with the Rouse matrix R. The overdamped dynamics of the beads follow the Langevin 
equations 

BU 

Ik = --^ : + y/2>ykBT£ i , (2) 

with £j being independent delta-correlated Gaussian white noise, damping coefficient 7 and 
thermal energy ksT. We consider the cases of a free chain as well as of a grafted chain with 
g (i) = 0. The Rouse matrix reads 



/ 



R 



(2-e) 
-1 



-1 







\ 



(3) 



2 -1 

v -1 iy 

with e = 1 for the free chain and e = for the grafted one. We pass to a dimensionless 
time, t = n/'yt, and neglect in the following the tilde in our notation. As long as the 
chain is intact, it is described by the standard Rouse model of polymer dynamics whose 
corresponding interaction strength reads k = 3k bT ' jb 2 (b 2 being the mean squared length of 
a single bond) [4,5]. 

The system of coupled monomers is sketched in Fig. [T] for a free chain. Superimposed is 
the equilibrium distribution of e«, 



^eg(ei) 



cxp 



20o 



, with 0o = (ej ) 



(4) 



which is the same for all i. 

The problem of thermally activated breakdown can be casted into a first passage problem 
of a reaction coordinate = q; L — q^i over a barrier of height AE = U (e^) — U (0) situated 
at Cb- In our harmonic model the reaction is assumed to be irreversible and to take place 
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Figure 1: One-dimensional chain of iV+1 monomers connected by harmonic springs. Superimposed 
is the interaction potential with a barrier of height AE at eb = qi-\ — qi. The activation rate is 
assumed to depend on the barrier height, the position of the bond in the chain and the system 
size. 

once the reaction coordinate reaches Cb in the harmonic potential well. Thus the absorbing 
boundary at eb introduces a cut-off of the harmonic potential. 

In order to study the thermolysis of the chain, all N bonds are assumed to have a cut-off at 
tb and the chain is broken as soon as the first reaches the barrier. For the systematic study 
of the mean first passage times (the inverse activation rates) T m f p (i) of individual bonds only 
one of them is breakable (the one under study, with a cut-off at e^), the remaining N — 1 
bonds are described by perfect harmonic springs. 

First, we consider numerically the thermolysis of the whole chain. Thus, the set of coupled 
equations ([2]) was integrated by use of a Heun integration scheme. Averages were performed 
over an ensemble of at least 10000 trajectories. Initial configurations were generated using 
the equilibrium distribution of the e^. 

The survival probability of a single bond in the chain is given by 



with the breakdown rate U{ equal to the activation rate of the bond i over the barrier, which 
on its turn is proportional to the inverse mean first passage time to e&, i.e., v{i) = l/r m / p (z). 
The survival probability of the whole chain (assuming uncooperative activation of the bonds) 




(5) 



5 




0.07 



0.065 



0.06 



Q 0.055 

Q. 



0.05 



0.045 



0.04 



2 4 6 



10 12 14 16 18 20 22 
bond i 



(a) 



(b) 



Figure 2: Panel (a): The fragmentation rate of the chain as a function of its length. Shown is 
v Nl v e- Panel (b): Probability distribution of activation for a chain with N = 22 bonds. The 
barrier height is AE/ksT = 5. 



is thus 



N 



W N (t) =Y[Wi(t) = exp[-u N t] 



(6) 



i=l 



with the fragmentation rate of the chain 



N 



U N 



£"(0- 



(7) 



i=l 



For a set of N bonds with equal activation rates v{i) = v e we have v^jv e = N. Due to the 
coupled dynamics this scaling with N is shown not to hold true. In Fig. [2] we depict the 
numerically obtained activation rates for free chains of different length. The scaling differs 
drastically from the linear one for small chains and approaches asymptotically a slope of 
almost one in the limit of long chains. The rate is always below its value for the case of 
identically activated bonds. We conclude, that in longer chains especially the inner bonds 
have lower activation rates, or — in turn — larger mean first passage times. Hence a chain is 
expected to be activated with higher probability close to its ends. This is shown in panel (b) 
of Fig. [2] where we present the probability density distribution of activation as a function of 
the location of the bond in the chain. 

The advantage of this simple model is that as long as the chain is intact, it is described 
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by the standard Rouse model of polymer dynamics, and therefore represents sufficiently well 
what happens in melts and concentrated solutions. 



A. Normal modes of the free chain 

The set of equations (j2J) can be decoupled by transformation to normal coordinates 
For the free chain the normalized eigenvectors (normal modes) of ([3]) are 



^« = l/^cos((, + i)^ T ) , (8) 



and the corresponding eigenvalues can be found as 



^=4sin(^-A-) . (9) 



2 N+ 1 

They are inverse proportional to the relaxation times of single modes. The decoupled equa- 
tions of motions describe independent Ornstein-Uhlenbeck processes with relaxation times 

rl = i (10) 

for k = 1,...,N. 

Representing the monomers motion with respect to the normal modes, the monomers 
coordinate is given by 

qi (t)=xo + 2Y^Xk(t)cos(^ + ^j^-^j , (11) 
with {^^(t)} being the new coordinates 

^ t ) = FTT^ 9i(f)cos (( i + S^r)' (12) 

Due to the choice of normalization the dynamics of Xq represent the time-evolution of the 
chain's center of mass. 



B. Normal modes of the grafted chain 



For the grafted chain the normalized eigenvectors of (j3J) are 



2 / 2k- 1 , 

xdi) = - sin iir—-. — - , (13) 

W V2iV + l V 2JV + V 



and the corresponding eigenvalues can be found as 



7i 2k - 1 x 2 



A£ = 4sin . (14) 

k \2 2N + l) K ' 

Representing the monomers motion with respect to the normal modes, the monomers coor- 
dinate is given by 

2 ,\ ■ ( ■ 2k ~ 1 \ 
qAt) = ■. > x k (t) sin itx— — . 15) 



k=l 

The relaxation times of single modes are 

r 9 k = i (16) 

for k = 1, N. Note that longest relaxation time of the grafted chain,rf , is approximately 
four times larger than t{ . 



II. THERMOLYSIS AS DIFFUSION-CONTROLLED REACTION 



In order to calculate the mean first passage times we briefly recall an approach 



pu t 



forward by de Gennes [16(, based on the pioneering work of Wilemski and Fixman |14j . [15]. 
The probability distribution functions of the distances ei are defined as ^(e^t). Their 
dynamics follow the generalized reaction-diffusion equation: 

&P = -QiP, (17) 

with £ being the diffusion operator in the absence of reaction and Q is the sink operator 
describing the reaction. Here it is assumed that the presence of reaction does not affect the 
distribution of all other variable ej, j ^ i. Choosing the delta-function sink (the Smolu- 
chowski sink), i.e., Q(ej) = K5(ei — e^), we have in the limit of infinite sink strength, i.e., 
K — > oo, an absorbing boundary at e^. The aim is now to derive an expression for the mean 
first passage time over a barrier corresponding to the energy growth towards e& into this 
infinitely deep and steep "adhesion well" . Eq. fflTl) can be formally solved using the Green's 
function method. The propagator Q follows from the solution of the equation 

Cg( ei , e°; t-t )= 5( ei - e°)5(t - to) . (18) 

Q is the conditional probability of finding the bond with elongation at time t provided 
that it was at e° at to- The formal solution of Eq. ffl7|) reads 



^(e i ,*) = Ve,(ei)- / dt I de°g(ei,e?;t-t Q )Q(e?)il>(e9,t ). (19) 
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The fraction of bonds that have not crossed the barrier at e& is given by p(t) and obeys the 
following relation 

= / deQie^ie^t) = K^(e b ,t) 

3 ft (20) 

= Kip eq (e b ) - K 2 dt g(e b , e b ; t - t )ip(e b , t ) 



dt 



Q(t) = Q(e b , e b ; t) is called the memory function. Taking the Laplace transform of Eq. fl20|) 
together with p(0) = 1 and K — > oo, we have 



1 - sp(s) 



sG(s) 



(21) 



For times exceeding the longest relaxation time the memory function approaches the steady 
state distribution ip e (e b ). Thus, one usually introduces a function which vanishes when 



t — > oo 



Eq. (l2"Tj) then reads 



m 

tpeq{e b ) 

h(s) 



1 . 



(22) 



(23) 



1 + 

It was shown [6|, H0(, that the long-time rate constant s* can be obtained by finding the 
pole of the previous equation. Furthermore in the limit of large characteristic time scales 
(potential barrier height much larger than ksT), s* is close to zero and the mean first passage 
time r m f p is approximately given by 

1 



T, 



mfp 



Mo). 



(24) 



Since the propagator is related to the two point joint probability distribution of via 

Q(ei, e°; t) = ipi(ei, e°; t) / i[) eq (e^) , the mean first passage time becomes 

'ipi(e b ,e b ;t) 



7~mfp (^) 



1 dt. 



(25) 



Following 6| it is shown in App. [A] that for any one-dimensional harmonic chain the two 
point joint probability distribution is given by 



^i(et,e?;t) 



exp 



e? + (e°) 2 -2c t (t)e^ 
20 O (1-Q(t)2) 



(26) 



2vr0 ov /l -Ci(tf 

with the normalized autocorrelation function Ci(t) = (ej(t)ej(O)) / (e 2 ) which we will specify 
in the next section. 
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III. RESULTS 



A. 



Free chain 



In the free chain the time-correlation function of = g« — is 



2 k B T 
N+l K 



N 



km 



{ei(t) ei (0)) 



W sin 



iV+ 1 



(27) 



where we have used the fact that different modes are orthogonal. The normalized correlation 
function Cj(i) is 



with r k taken from Eq. (fTUj) . For N ^ i the sum can be replaced by an integral, furthermore 
N + 1 ~ N. Substituting k by I = kn/N we have 



Thus for N ^> i, i.e., for bonds close to the chain's end, the correlation function Cj(t) is 
no function of the system size. The integral in Eq. (|29l) can be solved in terms of error 
functions. 

Since the mean first passage time is a functional of the correlation function, which by 
itself depends on the location of the bond in the chain, strong differences in the scaling of 
the correlation function might cause drastic differences of the activation rates (times) for 
bonds with different localization. Thus we first present in Fig. [3j panel (a), the normalized 
autocorrelation function of a free chain of 99 bonds. Depicted are the temporal correlation 
functions for a bond at one of the terminals of the chain (i = 1, dashed line) and for a bond at 
the center of the chain (i = 50, dashed-dotted line). Superimposed is the correlation function 
of a single bond, i.e., in a dimer (solid line). At short times (t <C r{) all curves coincide, 
while for longer times the correlations decrease much slower at the chain's center compared 
to the terminal where the correlation time (defined in the sense, that the correlation function 
is markedly different from zero) is also larger than for the dimer. Thus the correlation time is 
increased by orders of magnitude for bonds at the chain's center and the typical timescale of 
relaxation can reach and even overcome the mean activation timescale. Hence the dynamics 
become strongly non-Markovian. 




(28) 




(29) 
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Figure 3: Normalized autocorrelation function Ci(t). Panel (a): Free chain with 99 bonds. The 
correlation function is shown for i = 1 (dashed line, terminal of the chain) and i = 50 (dashed- 
dotted line, center of the chain). Superimposed is the correlation function of a dimer, i.e., a chain 
consisting of only one bond. Panel (b): Grafted chain with 100 bonds. The correlation function is 
shown for i = 1 (dashed-dotted line, grafted terminal of the chain) and i = 100 (dashed line, free 
end of the chain). Superimposed is c\{t) for a single grafted monomer. 

Since the analytical expression of the mean first passage time given by Eq. ( 1251) is exact 
only in the Markovian limit, the theory is expected to work more accurate in predicting the 
times close to the chain's ends. However, as shown in Fig. HI it offers a qualitative picture 
that can explain observations in numerical simulations (i.e., it works qualitatively well in the 
whole range of i, which is typical for the Wilemski-Fixman approximation). In panel (a) we 
present the first passage times derived from Eq. (1251) for three chain lengths and a barrier 
height of AE = bksT. In panel (b) we show the outcome of Brownian dynamics simulations 
for the same set of parameter values. Qualitatively the outcome of the numerical simulations 
agrees very well with the theoretical prediction. For bonds at the ends of the chain the first 
passage times are smaller compared to the activation times of the inner bonds. For these 
bonds the theory is also quantitatively in good agreement with the numerical simulations. 
As expected the agreement becomes worse with enlarging distance from the terminals of the 
chain. However, both theory and simulations predict an increase of the activation time with 
increasing system size, even for the bonds located at i = 1 and % = N. The observed effect 
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Figure 4: Mean first passage time as a function of the bond position in the free chain. Panel (a): 
Mean first passage times obtained from Eq. (|25p with the correlation function given in Eq. (I28p . 
Panel (b): Numerically obtained first passage times. The barrier height is AE = 5ksT. 

is large already in relatively small systems (see panel (b) for N = 21: r m j p (10) is about 40% 
larger than T m / P (l)) and becomes even larger in longer chains. Thus our study of the mean 
first passage times of individual bonds reveals that the activation times are smaller towards 
the chain ends what in turn cause there a higher probability of fragmentation. 

The barrier height defines the intrinsic time scale of activation. In Fig. [5] we depict the 
numerically obtained mean first passage times over barriers of different heights and a fixed 
chain length. The higher is the potential barrier the weaker is the increase of the times for the 
bonds at the center of the chain. This illustrates that the observed effect of the dependence 
of the dissociation time on the bond location is of highly non-Markovian nature. For lower 
barriers the activation times are comparable with the timescale of correlations in the chain 
and the dynamics is non-Markovian. We can conclude that the growth of activation times 
with enlarging distance from the chain's ends is the larger, the longer the chain and the 
lower the activation barrier is. 
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Figure 5: Mean first passage times as a function of the position of the bond in the free chain for 
different values of the barrier height as given in the legend. The times are given in multiples of the 
first passage time in a single bond system (dimer). The chain has N = 5 bonds. 

B. The grafted chain 

In the free chain we observed an increase of the activation times towards the central 
bond of the system. In the grafted chain there is also one free end, but another end is fixed. 
It was mentioned in the beginning that the longest relaxation time in the grafted chain 
(corresponding to the first normal mode) is four times longer compared to this time in the 
free chain. Thus we may expect that the non-Markovian aspect (induced by long correlation 
times) plays an even more important role in the barrier crossing dynamics. 

In the grafted chain the normalized correlation of the i-th bond elongation is 

«<«> = (^TT)E r ' (- (*^t) " *■ " ^M}) 2 ■ (») 

with r| taken from Eq. (fl6|) . 

In Fig. [3], panel (b), we show the normalized autocorrelation function as given in Eq. fl30|) 
of a grafted chain of 100 bonds. We depict the temporal evolution of Cj(i) for i — 1, i.e., 
the bond at the grafted end (dashed-dotted line), and for i = 100, i.e., the bond at the free 
terminal (dashed line). Superimposed is the correlation function for a single grafted bond. 
At short times the correlations drop down earlier compared to the single bond situation and 
the bond at the grafted end. Note that the dashed line of panel (b) coincides with the solid 
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Figure 6: Mean first passage time as a function of the bond position in the grafted chain. Panel (a): 
Mean first passage times obtained from Eq. (|25p with the correlation function given in Eq. (I30p . 
Panel (b): Numerically obtained first passage times. The barrier height is AE = 5ksT. 



line in panel (a). Thus we conclude that for short times the dynamics of the bond at the 
free end resembles the dynamics of a dimer. We can infer, that for short times the dynamics 
of a bond at the chain's free terminal is the same, no matter whether the chain is fixed or 
not. 

In Fig. El panel (a), we present the first passage times derived from Eq. (I25j) for three 
chain lengths and a barrier height of AE = 5k bT. In panel (b) we show the outcome of 
Brownian dynamics simulations for the same set of parameter values. As in the free chain, 
for bonds at the loose terminal the first passage times are smaller compared to the activation 
times of the inner bonds. For these bonds the theory is also quantitatively in good agreement 
with the numerical simulations. The first passage times gradually increase with enlarged 
distance from the free end of the chain and grow substantially at the grafted terminal. The 
effect is overestimated by the theory. However, theory and numerical simulation are in good 
qualitative agreement. 

As for the free chain we present in Fig. [7] the numerically obtained mean first passage 
times over barriers of different height and a fixed chain length. The lower is the potential 
barrier, the stronger is the increase of the dissociation time along the chain. In the limit 
AE 3> k B T T m f p (l) — > r^f p while r m f p (i ^ 1) — > T^f p /2 (the value expected for a free 
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Figure 7: Mean first passage times as a function of the position of the bond in the grafted chain 
for different values of the barrier height as given in the legend. The times are given in multiples of 
the first passage time in a single bond system. The chain consists of N = 5 bonds. 

dimer) . 



IV. THE FREE 3D ROUSE CHAIN 



Let us complete our study and turn to the three-dimensional harmonic chain. A^ + 1 
beads are connected by A" harmonic springs. The ends are free and the chain's center of 
mass diffuses with an effective friction ~ N. In the three-dimensional chain system the 
two-point joint probability distribution of is given by 



?Me;,e°;t) 



2tt0 o 



1 



exp 



1 e? + (e' 



0\2 



2ci(t)ei ■ e\ 



20o 



i - Ci(ty 



(31) 



Q(t) 2 ) 3/2 

The derivation follows the same steps as shown in App. [A]and is given for example in 
After averaging over angles the distribution reads 

Q(t)e;e° 



2e 4 e° 



sinh 



Ci{tf 



exp 



Together with Eq. (1251) the mean first passage time reads 



T, 



mfp\ 



efc{t)y/l-Ci(ty 



sinh 



<k(t)el 



^o(l-Q(t) 2 ) 



exp 



ej + (e°) 2 
20o (1 - Ci(t)< 



(l-Q(t) 2 ) 



(32) 



I I dt . 

(33) 
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Figure 8: Mean first passage times as a function of the position of the bond in the chain. The 
barrier height is AE/ksT = 10 . The chain length is N = 9. 

with d(t) given in Eq. ( 1281) . 

In Fig. [8] we compare the mean first passage times of |ej| derived from Eq. ( 1331) and 
from Brownian dynamics simulations for a barrier height of AE/ksT = 10. As in the 
one-dimensional system, there is a qualitative agreement between theory and numerical 
simulations. Even for a relatively short chain consisting of N = 9 bonds, there is an increase 
of the activation time of about 8%, which will be even larger for longer chains and/or lower 
activation barriers which are of relevance in biological systems. 

V. SUMMARY 

Let us summarize our findings. We studied the thermally activated fragmentation of 
a homopolymer chain. It was shown that the fragmentation rate of the chain follows a 
nonlinear scaling as a function of the number of breakable bonds in the system. Bond 
breakage happens with higher probability at free chain ends. Studying the activation times 
of individual bonds we focused on the impact of their location in the chain and the length 
of the latter. It was found that towards the center of the free chain as well as towards 
the grafted terminal of the fixed chain the activation times increase substantially. The 
theoretically predicted effect is qualitatively confirmed by Brownian dynamics simulations 
both in one- and three-dimensional systems. The observed effects are large and therefore 



-■-numerical simulation 
— Eq.(33) 
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the framework presented in this article may help to interpret real experimental data. 
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Appendix A: TWO-POINT JOINT PROBABILITY DISTRIBUTION 

The two-point joint probability distribution is given by 
V>»(e<>e£;*-*o) = J - J dqdq S(e l -(q i -q i - 1 ))5(e° i - (g° -g,°_i))^o (q, Qo; * ~ *o) • (Al) 

Since ipo — the joint probability distribution of all of the chains coordinates — is a multivari- 
ate Gaussian, the integral in Eq. flAll) yields a multivariate Gaussian {20]. The general 
representation is 



ipi(e h e1;t) = N exp 
with iV being a normalization prefactor. From 



\mi{t)el - \m 2 {t){ey - ~m 3 (*)e ie ° 



(A2) 



de°V>i(ei, e"; i) = ^ eq (t^) 







(A3) 



(A4) 



and the normalization restraint of the probability density distribution we obtain mi = m 2 
m and N = m? — m|/4/(27r). We derive 



1peq{ei) 



III 



4m 



2tt 



exp 



m — 



Am I 2 



(A5) 



Furthermore we derive (e 2 ) = <fio = l/(m — m\j '(Am)). From 

■bt 1 



Ci(t) 



(A6) 



it follows that = — m 3 /(2m) and subsequently 1/m = </> (l — Cj(t) 2 ). Finally we have 



1 



27T0 OV /1- Cl (t) 2 
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exp 



e 2 + (e°) 2 -2c t (t) e ^ 
20 O (1-Q(t) 2 ) 



(A7) 
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